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By H. Latzko 



The theory of heat . transfer from a solid "body to a liq- 
uid stream could he pyeisented previously** only with limit- 
ing assumptions ahout the movement of the fluid ' (potential 
flow, liaminar frictionai flow)* (See references 1, S, and 
3.) 

Por turhulent flow, the most important practical case, 
the previous theoretical considerations did not go teyond 
dimensionless formulas and certain conclusions as to the 
analogy hetv/een the friction factor and the unit thermal 
conductance, (See references 4, 5, 6, and 7,) In order to 
ohtain numerical results, an experimental treatment of the 
pro'olem was resorted to, which gave rise to numerous inved^- 
tigations hecause of the importance of this proTjlem in many 
"branches of technology. Hoxvever, the results of these in- 
vestigations frequently deviate from one another. The ex- 
perimental results are especially dependent upon the over- 
all dimensions and the specific proportions of the equipment. 

In the present work, the attempt will he made to devel- 
op systematically the theory of the heat transfer and of the 
dependence of the unit thermal conductance upon shape and 
dimensions, using as a hasis the velocity distrihution for 
turhulent flow set up hy Prandtl and "Von Karm^n, 

*"Der Warmetltergang an einen turhulfenteh Pltlssigkeits- 
Oder Gasstrom,'' (Ahstract of a Dissertation presented to the 
Phil. ' Faculty of the Univ, of Vienna.) Z.f.a.M.M. , vol. 1, 
no'. 4, Aug, 1931, pp^ 368-290. 

♦"'As long as the velocities remain much "below the veloc- 
ity of sou?xd, comprassihle fluids (gases) and incompressible 
fluids follow, as is known, approximately the same laws of . 
flowT"'ffi'&t'W6"i''ei ±n the f ollOAif4ng.,v, 4hi&,4!xpression for the 
flow of fluids v/ill be used for actual liquids'as well "as 
for gabes. 

HOTS; Translation received from Univ, of California, 
Berkeley 4, Calif, 
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2. 



1. HYDilODYirAin 0 PRIITOIPIiBS 



First of all, the resTiltg of the Prandtl-Karman theory 

(reference 8) which will "be used constantly, v/ill he set 
forth, 

For the distrihution of the shear stress in the immedi- 
ate -vicinity of the wall, dimensional considerations (see 
note 1 in the appendix) yield the expression;. 

T = 5- To p n y ^ .. (2.) 

where the symhols are: 

u velocity in the direction of flow 

. shear stress at the wall 

|j. ahsolute, viscosity 

V = Ia/Pi kinematic viscosity 

B constant (note 2), (8,82)* 

y distance from wall 

p mass density 

In the same region, when the shear stress at the wall 
is assumed as known, the velocity follows from the equation 
(note 3): a,j \, . 

/"^o\ 7 /y\ 7 

There are two methods (note 4) of ohtaining the distri- 
hution of velocity and shear stress for the entire region of 
the fluid. Either one starts from equation (2) and sets 

u (y ) = y , (^Ao + Aiy +. Agy^ + . . . ) ( 2a) 



*The value of the constant corresponds to the equation 
for the velocity . distrihuti on , virhich is used helovr. Compare 
eq.uatio.n (8), 
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in v;hich the consts^nt Aq is determined from the require- 
ment -t'hat equation (3a) is transformed ,irut.o equation , (3) at 
small values of y, or the basic equation (l) for the shear 
stress transmitted Isetween the individual layers ban "be ex- 
tended and equation (la) can oe written 

r 

. T K I^Y(y) 

in which Y leust change to y in the vicinity of the wall. 

With this TDasic equation the velocity field for turTsulent 
floxiT can he calculated as long as no separation from the 
Tsoundary walls takes place, 

?or the special case of flow in a right circLilar cylin- 
der it v/as shown hy Von Karnan that the experimental results 
on tlie velocity di stri"buti on can "be reproduced with siiffi- 
cient accuracy if the function Y (y) is made, (called the 
influence function). 




where y = r - y, 

For time and volume invariant, the velocity distribu- 
tion is then: 



bu 



(la) 




and finally, the relation of the maximum velocity at the 
axis of the tuhe Umax ^^"^^ .average velocity v in the 

cross section is 

= S' ^inax (4) 



2. TUilBUlBITS THSHiiAL COITVEGJIVITY 



In the follov/ing, the transmission of heat "by matter 
only is considered, conseqixently limiting the study to a 
temperature region in- which the amount of heat carried ofjf 
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"by radiation ie negligilJle in comparison with that carried 
of f Tjy'^^rticXe'?" of yurthermo-re-, the vel oci ty - dis- 

triTJution shall Tie affected only "by external conditions; 
that is, the influence of the temperature field upon the ve- 
locity field is disregarded. At relatively great velocities 
of flow vrhere. the notion is tur'buleht, the resulting error 
iaeed not "be taken into account as long as the differences in 
density in the cross section, caused "by temperature changes, 
are not too great. 

Corresponding to the ideas taken from those on the con- 
duction of heat in solid "bodies, a distinction is likely to 
l3e. made in the case of heat transfer in fluids, in general 
between the thermal o ondtxctivity , which descrilses the heat 
transported Tjy molecular movement, and the so-called thermal 
convection - that is, transfer of heat by movements of the 
mass. Ihe order of magnitude of the carrier of heat is thus 
used as the basis for distinction, A somewhat different 
mode of consideration, which pushes into the foreground the 
nature of the motion of the carrier of heat, seems, however,' 
to he more advantageous "both for mathematical treatment and 
for comprehension of the process. Accordingly, "by "thermal 
conduction" in liquids is understood the transmission of heat 
"by the random motion of the molecules, as has "been repi*e- 
sented "by the concepts of the kinetic theory of gases. It 
then will "be regarded as characteristic of the molecular 
movement of heat that it is a pure function of temperature 
at a fixed pressure and a fixed density of the fluid and 
especially that it is not dependent upon the state of motion 
of the fluid. Thermal convection, on the other hand, 
shall signify the transfer of heat which results v;hen the 
motion of the particles is directed. In many textbooks of 
physics free convection is considered as the origination of 
a natural flpw produced by differences of density under the 
inflxxence of the force of gravity, IDhis concept,- then, is 
contained in the preceding definition. 

In the case of laminar flow, all the heat transfer can 
be accounted for by the foregoing concepts, For turbulent 
flowj however, one manner of heat transfer is still unmen- 
tioned. As i s known, s-t eady- stat e turbulent flow is repre- 
,s,e.nt^e4 a.p. having .(^.fe .eaoh,pP.int a oertaiji average velocity 
vector upon which is superposed ano'the'r velocity vector, 
varying in . directi on. and magnitude, having an average value 
over a sufficient span of time equal to zero. According to 
Ton Karman this kinematic picture can be described more ex- 
actly by the representation that vortox filaments with a 
random motion float in ths bulk of the fluid, which moves 
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along with a fixed time-average velocity distriTiuti on, 0?he 
"^d^efflgatS'' of the vor-t-ex filament^ , , as well those of the 
molecules, obey laws of statistics, The fiuct-ixating veloc- 
ity vector (time average is zero) is then defined at a point 
of the fluid "by the circulation and "by the relative position 
of all the vortex filaments, 

This concept leads to introduction, apart from the usu- 
al thermal conductivity, which appears as an expression for 
the statistical law of molecular motion, of a conductivity, 
of turbulent motion which expresses the statistical influ- 
ence of vortex motion upon the transfer of heat. It then 
will depend primarily upon the state of motion of tho fluid, 
which is especially influenced "by the nature of the boundary 
surfaces. 

The method of accounting for this phenomenon by intro- 
duction of an increased conductivity for turbulent motion is 
known. Several authors have proposed different basic equa^ 
tions in which the increased conductivity is regarded as an 
empirical function of the velocity. Recognition of the true 
circumstances was partially clarified by the considerations 
of Reynolds and Prandtl, Both began with the idea that tur- 
bulent friction and turbulent heat transfer are analogous 
proces^ses, and that the same mechanism which in the first 
case causes a "momentum transport" leads to transfer of heat 
in the second case, Reynolds (reference 4) in an intuitive 
manner, according to this consideration, went directly from 
the friction factor to the unit thermal conductance in cir- 
cular tubes and compared, as it v/ere, the integral procosses. 
On the other hand« Prandtl (reference 7) sets up the exact 
conditions under which a directlrf analogous conclusion is 
permissible; he shows that in certain cases the temperature 
field is an exact image of tho velocity field, so that 
knowledge of the motion permits direct conclusions about tho 
thermal field. However, ho shows that this is clearly not 
the case for right circular tubes, so that co.ncl-asi ons can 
only be drawn as to the form of tho relation between tho 
different parameters, since numerical results cannot be ob- 
tained. Recent advances (sec. l) in the mathematical repre- 
sentation of turbulent flow and the velocity distribution 
corresponding to it now make possible a more exact expres- 
sion of the "elem'entarj' law" for turbulent heat exchange, so 
that the following statements start out from Prandtl 's re- 
sults in two directions, by first of all furnishing numeri- 
cal results and then allowing a mathematical consideration 
of the different arrangements in v;hich there exists no spa- 
tial constancy of the velocity and temperature fields. This 



HJlGA. m ITo, 1068 



6 



malces possiTjlis a detailed dlactission, of .experimental results 
to explain individual deviations* 

3. PUNDAMBirTAL LAW OP TURBULENT HEAT BXOHAN&E 



First of all, the processes of thought applied to the ^ 
kinetic theorjr of gaaes, which lead to the differential law 
of internal friction and heat transfer "by random motion of 
molecules, will Tie applied to the case of turhnlent exchango. 

Consider a layer at a distance y from the wall; there 
the average velocity u prevails in the direction of flow 

and let; |wi he the average absolute amount of the velocity 
perpendicular to direction of flow. Then the average velo- 
city of flovr in two layers at a distance *~ from the layer 

3 

y under consideration (where x is a kind of "mean path") 
i s : 

u ± ^ 5. 

by 2 

The momentum transport per unit of surface perpendicu- 
lar to the average flow is given, introducing a proportion- 
ality factor P which depends upon the nature of the coher- 
ent parts of the fluid and the formation of the mean value 
with respect to time, ty the expression: 

Ppw X = T (5) 



and is equal to the shear stress T at y» 

If 0 is the heat capacity of a unit volume, then, on 
the other hand, the heat transport q. per unit of surface, 
likewise perpendicular to the average flow, is given hy : 

^ = pOw M X (6) 

This states that the same fluid particles which produce 
the shear stress T hy their transmission of momentum, also 
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transmit the heat ^ l?he proportion- of .hea,|,.i;rai:i_sinitted can 
"be calculated "by a kind of counting of these particles* 

One such "counting" is given, as is 'easily seen, "by the 
product Pwx (called the coefficient of turbulence). The 
cboff icient of turhulenco, together with th© constant 0, 
represents an expression for the statistical law of heat 
transfer in the case of turbulent tlov just as does the con- 
ductivity X in the case of no flo\7» 

She coefficient of turbulence can be calculated from 
previously obtained knowledge of the state of flow, 

From equations (l) and (la) 



T = pcp(y) 1^ 



from which (note 5) 



7 T, 



3/, 



Pwx = cp(y) = ~ y% 



(7) 



The basic equation (l) expresses the total effect of 
the molecular conduction of jnomentum (internal friction) and 
of the momentum transport by eddy convection. Correspond- 
ingly, the basic equation itself, as well as the velocity 
distribution originating from it, is to be regarded only as 
an expression which becomes asymptotic at very great 
Reynolds numbers, where the effect of the molecular conduo- 
tion of momentum is small in comparison with the second part 
of the friction mechanism - that is, eddy convection. How- 
ever, it has been found that the proportionality between t 

and V * is a very good approximation even at values of the 
Reynolds number which correspond to about five times the 
critical velocity. From this, it is concluded that the sta- 
tistical laws for the molecular and eddy transport of momen- 
tum can be represented. to a good approximation, even at 
moderate Reynolds numbers, by th© general expression (l)^ 

By referring to equation (l) for caiculation of pwx, 

it is assumed that all the heat transfer can be expressed 
also by a general statistical law which summarizes molecular 
and eddy processes. It has been assumed, therefore, that 
there exists in the molecular processes the same proportion- 
ality between momentum and energy transfer as exists in the 
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eidy prJ3J?A8Sj8$;, that is, it is assumed that the ratio ■be- 
tween ^ and |x is the saine as that T5ett?een= 0 and p. 

The error committed is negligilsle for gases, as is 
shov/n "by the following consideration. On the one hand, the 
portion of heat carried over hy pure t^irhulent convection is 
several times that transferred "by molecular condition, as is 
sho.v/a "by a comparison of unit thermal conductance for lami- 
nar and turljulent flow; on the other hand, the ratio ^ 

lies "between 1^26 and 0.97 according to the number of atoms 
in the gasj that ivS, the molecular mechanism of condition of 
momentum (internal friction) and that of conduction of heat 
are essentially similar, Hence, for gases and superheated 
steam, practical and quantitatively correct results can he 
expected, from the calculation. 5he follox^ing derivations 

are to Tjo understood in this sense. The case where — ^ 

Cp. 

differs greatly from unit-y will he referred to once again at 
the conclusion of the vork* 

By consideration of equations (6) and (v), there is oh- 
tained for the total amount of heat q transferred throufjh 
a unit surface of a layer at a distance y: 



7 /T 



In practice, the limiting value of q^ for y = 0 - 

that is, the amount of heat going out of the wall per unit 
of surface w will he calculated as follows: The velocity u 
is represented hys 

>i(y) = y^'^7 I^Ao + Ajy + Agy^ + , , .j 

The shear stress ^ has a fixed limiting value for 
y - 0, and is a regular function of y in the vicinity of 
y » 0, Therefore, T can he developed as a pov/er series in 

y. 
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Inserting , 

- — cty . 

^ 7 -^0 + ^IJT + . . . 

then developing this fractional expression according to 
powers of y in the region y = 0; 

finally, considering eq.uation (S), yields: 

p ^0 y=o lay J 

4. HBAO? SZCHAiTCJE liT TUBES 



When there is steady-state flow through a tu"be, two 
regions can he distinguished: 

1, Fully developed flow state - that is, one in which 

the velocity-profile remains similar along the direction of 
flow 



2, The hydrodynamic calming length at the entrance to 
the tuhe 

Assume, for example, that the fluid flows into the tuhe 
through a smooth passage from a large reservoir; then at the 
inlet cross section the streamlines will have approximately 
equal velocity. On progressing further, the layers near the 
wall v/ill he retarded "by friction until the constant (with 
length) velocity profile, iirhich corresponds to the steady 
state, has "been developed, Shis part of the tuhe is often 
called the entrance ' secti on. 

In the following sections the temperature field and the 
heat transfer in the tuhe are calculated for the case where 
a, temperature distriljution for the entrance section is given 
heforehand and the v/all temperature is kept constant along 
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the direction of flow. Separate solutions are set up for 
the tw^' "r^egions TtteTTt'loned, "but hy a continuous transition 
from the first solution to the second, they can satisfy'the 
general function through summation of the partial solutions. 

5. KEAT TEANSPES PGR IHE CASE IN WHICH THE VELOOIIY 
DISTRIBUTIOK HAS ?BEIT ESIABLISHBD AT THE 
BWTHANOB TO THE THEBMAI SECTION 



In order to set up the differential eq_uation for the 
temperature field, an element of volume, "bounded on the 
sides "by two concentric cylindrical surfaces, parallel to 
the walls of the tube, and TDOunded on the "ends "by the cross 
sections perpendicular to them, is considered. In order to 
complete the representation it is assumed that a warm fluid 
flows through a colder tu"be; that is, the flow of heat shall 
"be from the fluid to the wall, Furthermore, the constant 
temperature of the tu"be is set equal to 0, so that the fluid 
temperature is the excess temperature above that of the V7all» 
However, since no assumption is made which distinguishes one 
direction of heat flow from the other, all relations are 
valid when ^ changes its sign. 

If z =! coordinate of the direction of flow 



from equations (s) and (4), u is replaced hy (note 7) 



y = distance from the axis (note 6) 



C 3 heat capacity per unit of volume 



then the heat "balance for the steady state gives: 




(10) 




while q follows from equation (8); then, also considering 
(la), (3), and (4) 
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• ■ - ■ ■■ - ■ ■ . V V 0 

^ 0.19-9 

(2r) 

« 

If this epilation (note 8) is introduced into epilation 
(10), there is finally o^btained as the differential equation 
for heat transfer when the conditions of flow are hydrody- 
namically complete: 



8 ^'UZv)"'-r (lla) 



where 

K 



^ 0,199 V^^"^ 



The "boundary conditions are; 
I, d » 0 for y = r. 

H, hA 8 0 for y = 0, "because of the universal symmetry 
by 

III, The radial temperature di strihuti on must "be given 
for z a 0. 

Since the fluid temperature approaches asymptotically 
the. temperature of the wall, as the tuhe-longth increases, 
then the solution is of the form 

^ - g(y)e-^'' 

If this expression is inserted into the equation, then 
there is ohtalned for the function q the ordinary differ- 
ential equation 

which, after elimination of the fractional exponents Ijy the 
transformation 
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2^ " 




1 - ^ J T = X 



.r 



'becomes (note 9) 



J 



where (note 10) 



(12) 



6/ 

w = 49kK 



(I) 

The Tjoundary conditions are now: 

I, e = 0 for X = 0 
de 

II, — is finite for x = 1 
dx 

An approximate solution is ol3tained liy means of the 
Ritz method (note 11), when the proljlem is changed to one in 
the Oalciilus of Variations (note 12); that is, as can Ise 
verified easily (note 13) : 

-1 s 

(l - x'^ ) ( " wx'^'g^l dx = rainimiim (iS) 
\dx/ J 

u 

with the "boundary conditions I and II as supplementary con- 
ditions. Here the problem is one of finding the "character- 
istic values," since equation (13) will have solutions which 
also satisfy the houndary conditions, only for fixed values 
of w, Suhstituting for g! ' 

. g(x) = giPi(x) + g3P3(x) + g3p5(x) + . . . (14) 




in which gj^ , gg, are undetermined coefficiente and 

-^3 » and so forth are the Legendre spherical functions 
(note 14) of the first kind, and taking only three terms 
first, results, for w at the minimal conditions, in an 
eq.uation of iTxe third degree the roots of which are 

Wi = 8.712, Wg =; 164.36, W3 a 1700.40 (l4a) 

The characteristic functions are normalized, in con- 
trast to the customary procedure, so that 
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(15) 



in order that the temperature at the axis of the tuhe will 
Tiecome unity; e(y) is, therefore, the ratio of the tempera- 
ture at the point under consideration to the temperature at 
the axis of the tnlje-. Thus the first characteristic func- 
tions are: 



^III " 2.6552 Pi - 6,1589 P3 + 4,5037 P5 

The choice of development according to spherical func- 
tions must he justified. Since it is a minimal prohlem, the 
exact characteristic values can only he less than the ap- 
proximate values. The magnitudes of the characteristic val- 
ues which are obtained according to the choice of the hasic 
series equations for the function to he varied, form a suit- 
able criterion for the validity of the approximation. How, 
it can "be seen that, compared to a simple power equation in 
X, as well as several Courier developments, the "basic equa- 
tion in spherical functions leads to the least characteristic 
values. As to their behavior on making further approxima- 
tions, the first three appro;cimati ons for the first charac- 
teristic functions furnish , successively , for example, the 
values 8,75, B,67, 8.71; consequently the convergence of the 
procedure ought to he satisfactory. 

In order to obtain equally good results for the other 
characteristic values, further approximations must naturally 
be made; the third characteristic value, in particular, will 
agree only in magnitude in the case of a three-member basic 
equation. As will be seen, however, this has only a slight 
influence upon the results, 

The particular merit of the spherical functions for the 
problem in hand also can be demonstrated by the following 
simple consideration, which can at the same time dispel 
doubt caused by the increase of the second coe.ff icients in. 
the second characteristic function,. Considering figure 1, 
it is 9een that the characteristic values themselves show 
great similarity to. spherical functions. If a form, like 
that represented by equation (16), is now set up according 
to functions which are identical with th^ characteristic 



0.9703 Pi + 0.0212 P3 + 0,0085 Pg 
0.7312 Pi + 0.9665 P3 + 0.7647 Pg 



(16) 
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values, then the coefficient the index of which is equal to 
the ordinal ' humher of'the characteristic value Is ecomes actual 
to unity,. .all others being equal to zero, The characteristic 
value concerned is already represented exactly by one member 
of the development,. If developed according to functions 
which are not identical with the characteristic values^, but 
which have, however, a certain ' similarity to them, then, in 
the development the coefficient having an index eiqual to the 
ordinal number will be slightly greater. Clearly, such a 
development will islosely approximate, with relatively few 
members,, the function to be represented. 

If the values of w from (l4a) are inserted into equa- 
tions (lla) and (l2a), then for the coefficients of the ex- 
ponents 

1 % 1 i % 

ki = 0.1510 ^ (:^) ; ks = z.m e[ (^) ks = zs.^ '-^ ^ (17) 

The complete solutibn of the partial differential equa- 
tion (13) can be \irritten as a development according to char- 
acteristic functions 

d ^ aigje"''^^ + asg^je'^^s^ + a^gjue'^^^ (l8) 

where the coefficients are to be determined so that the pre- 
scribed temperature distribution is fulfilled for z = 0, 
The calcxilation will be carried out first for the case of 
uniform temperature distribution at the initial cross sec- 
tion. Therefore aj ... ... &^ first must be determined so 

that d (y) must be as close to ^ = 1 as possible.. Every 
other temperature then follows with the aid of a multiplica- 
tive constant,.. The least square error, yields the values:- 

ai =5 1.129, as, =i -0.180, a^ = 0.048 

ft 

Thus the final equation of the temperature field for 
turbulent flow in a hydr odynamically complete state in tubes 
for the case- where. th.ery;,ni,,fO|;,m temperature prevails in 

the initial cross section, is 
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f V \ z 

... . . -.<^,.i§i(^ y -J 

§f§p<'ia§9 e, .,[0,9544 .X 5- 0 

«Q,180 e • ■ C-0.7472 x- 4.275 6.023 x**] 



>(19) 



•a 9*4 a' -"^ ' ^ 



+ 0.048 e .^[2.0.34 X- 54.80 + 35,47 x*] 



J,n v/liich the similar powers of x are collected from the P 
^Jjegendre spherical functions). 

It is also recognized now that an error in the third 
oJ;j^ara"cteristic value and in the third characteristic func- 
•y^on is of slight significance; even if the third exponent 
a^^(^uld "be still somevAat smaller, the third characteristic 
fuHijQttion dies out several centimeters from the "beginning of 
thQi^ thermal effect, the error having no influence upon the 
remaining part of the tuhe. A fourth approximation always 
can he calculated. 



6, DISCUSSION 01" RESULTS AID AGREEMEErT WITH EXPERIMBIIT 



By reference to figure 1, the temperature distrihution 
over the cross section (of the fluid stream) can he discussed. 
Por 2=0, a sauare distribution was assumed; that is., the 
fluid enters with a uniform temperature over the whole cross, 
section. 

In the interval "betv/een 0 and 0,8 for ~, the uni- 
form temperature is represented to a maximum error of ±2^ 
(per 1000) hy equation (l9), 

y 

Tor — "betxireen 0,9 and 1 -there is a sharp temperature 

decrease since only thr ee' terms 'i^Bir g ednsid'erigd, Similar 
situations also exist in reality, since the layers near the 
wall will undergo a change in temperature, due to radiation, 
"before making direct contact with the wail,. On moving far- 
ther alone the tu"be, the temperature, gradient at the wall 
levels out more' and more; the so-called final temperature 
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,di s tritviti on , which is represented "by the first characteris- 
tic function only, is reached when the sec6hd characteristic 
function has died away. From then on, all temperature pro- 
files remain similar since all temperatures decrease in the 
z^direction according to the same exponential function; the 
expression "final temperature distribution" is to te under- 
stood in the ahove sense. As shown "by eq.uation (l6), the 
first characteristic function and, consequently, the final 
temperature distritiution, differ Tjut little from the velo- 
city distriljution in the hydrodynamically complete state,, 

■ With the help of the known temperature field the point 
now is reached where all the questions ahout the heat trans- 
fer can "be ansv;ered. For example, to calculate the unit 
thermal convective conductance a, the ratio is set up of 
the amount of heat transferred per unit of wall surface for 
the mixed~mean temperature difference at the cross section; 
that is, 

m 

According to equation (9), q- is given hy 



(30) 



0..176 v^^^ Ov"^^' 



' I- [t^ ^^'] 



and the average temperature is defined "by the equation: 

r 

Hence, the expression (note 15) for ot. isf 

a=:C.0346 vC ('JlV'^ 1. 078 e-^^i^+0.1S4 e"^^^ 0. 98 0 e'^a ^3^^ 
Vvd/ 0,970 e~^i^+ 0, 024 e^^s^ + 0, 006 e'^s'^ 



The analogous result for laminar flow was calculated "by 
]^usselt, (See reference 6.) 



Figure 2, which is calculated for the special value 
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d \v.d/ 



% 

= 0, 037, sliows the va^xatibn of ' "i;he unit tHetaal 



convective conductance A^itli distance into the tuTDe,/ At 
2=0, a is infinite, then, corresponding to the decrease 
in the temperature difference at the wall; it decreases,, 
though considerahiy faster than it does when the flow is 
laminar, finally approaching a minimum value %iin» •^1" 
though this least value is independent of the velocity in 
the case of laminar flow, for turl>ulent flow equation (2l) 
is changed into the form: 

V . • 

ttj^ijj = 0.0384 vO (-^^ * (21a) 

Equation (21a) is analogous to the equation developed 
"by Seynolds, • Since, as mentioned already, the eddy heat 
transfer in turhulent flow exceeds the molecular one "by a 
multiple, it seems justifigihle that only those magnitudes 
which are determinative for the condition of flow and also 
for the eddy transport of heat should appear in the formula 
for the unit thermal convective conductance, -iPhese are mag- 
nitudes V, d, and y or C, . The variation with tempera- 
ture depends upon the values of the kinematic viscosity V, . 

If the relation for gases ^ 'sr 1 is considered, it is oh- 

OiJ. ~ 

served that eq_uation (31a) likewise agrees in form with the 
dimensionless formulas of Nusselt and Prandtl, • -(See rofor- 
ences 2 and 7, respectively.) 

Ihe existing experimental material is not sufficient,, 
unfortunately, for an exact test of these results, since 
average unit thermal convective (Sonductances were always 
measured and the "entrance sections" were not chosen long 
enough so that in the measuring length a hydrodynamically 
complete state with a temperature profile which remained 
similar could have "been attained with certainty, . In most 
cases the point at which the thermal effect hegan cannot 
even he determined. Ohvi;ously this is an indication that 
the experiB^en,ters possih|.y did not h|uve_, a clear concept of 
the influence of the arrangement upon the results of the 
measurement. 



Nusselt, in a short series of experiments*, ^connected a 
♦Research work puhlished in reference 7, tahl© 6. 
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piece of tuning 2 meters long in front of the actual experi- 
Sental geBtion, According to the calGi.ilati.ons-, tirhioh are 
discussed in the next section, the state of flow was cer- 
tainly complete. Besides, since this' entrance section and 
the heated experimental tu^be "both \irere made of Tsrass and 
were joined firmly to one another metallically, the added 
length of tu'bing likei/ise heated, at least in the part 

directly connecting with the experimental tu^be. 3?he first 
point of temperature measurement was, on the average, ahout 
15 centimeters downstream from the "beginning of the test 
section-, so that it can safely be assumed from the results 
of equation (Sl) that the unit thermal convective conduct- 
ance had reached its minimum value, Ihese experiments v/ere 
checked, using the equation for oi^irx- 5Ph-e results are pre- 
sented in the following tahle (p, air -pressure, V , unit 
\*eight density of air): 



Exper- 
iment 

num- 
her 




I 


Pm 


'Ym 


V 


*^eas 


<*calc 


Differ- 
ence 

(percent) 


95 


39.0 


0,6133 


1,161 


1.273 




19.29 


20.09 


U.lU 


96 


37. S 


.62U5 


1.167 


1.2S5 


5.75 


2l4-,95 


25.36 


l,6U 


97 


.3U.2 


.6368 


I.16U 


1.255 


S.29 


32.75 


32.57 


-.5 


9S 


31.5 


,6it3g 


ia63 


1.307 


13.06 


H6,S 




1*3 


99 


35..6 


1,0590 


1.16H 


1.291 

■ 


21.06 


65.3 


67.51 


3.35 


100 


32>1 


1,1300 


1.167 

■ 


1.309 

■ 


2k,Q3 


73.0 


75.25 


■ 5.0 



Experiment 95, at a Reynolds numljer of 6100 (about three to 
four times the critical velocity), was near the limit of the 
region of the validity of the ahove-developed theory, Con- 
sidering the limits of the aoc.u,i»,aoy of measurements of this 
kind» the agreement seems to be absolutely satisfactory. 

In the following sections, the heat transfer in the 
calpiing length of a tube ;i(ihat is. the heat transfer for the 
hydrodynajnically incomplete -state) .yfill be investigated. 
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Inasmuch as a solution for the velocity field in the calming 
length pf a tube has not been given previously, it must 
first of all be determined, . , - 



7» THE VELOCITY FIELD IN THE "CALMING LBNaTH" 



In order to obtain an approximate expression for the 
velocity field in the calming length of tubes under condi- 
tions of tuj^bulent flow (see reference 9, for the case of 
laminar flow), the momentum consideration introduced by 
Von Kdrman (reference 8, pp. 235 and 356) will be used. 

Consider a longitudinal section through the beginning 
of the tube. At A, fluid from a large reservoir flows into 
it with a uniform velocity distribution, The layers near 
the wall will bj3 retarded under the influence of the viscos- 
ity, and the thickness of the layer, in which the shear 
stress is transferred (shown by the shaded lines in tho fig- 
ure) will increase until the two boundary layers meet. From 
then on, with the insertion of a short transition region, 
the velocity distribution over the cross section will remain 
constant* 

Hence, it is assumed that there is at the beginning of 
the tube, a region in the interior of the flowing fluid 
where viscosity can be neglected. For this region, the va- 
lidity of the Euler equation, formulated for frictlonless 
flow, is assumed. 

If the steady state is assumed, a balance on an element 
of the boundary layer is considered, which has a ring-shape 
structure; abed in figure 3 represents a cross section 
of this element* 

Let 

Q a: the volume flowing through the cross section in the 
boundary layer per second 

J a the transport of momentum per second in the direction 
- of flow through th,e ^cjeg3Sjt$.%Q.;^%ojo.a.l surface 

U = the veLocity of free stream 

u = velocity in the boundary layer 

*L, Schiller, who most recently studied experimentally - 
a theoretical explanation also was published - the problem 
of the entrance section^ treats only the laminar caso. 
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y... 




average .ve.lOG.ity ovet the cross section 


8 


s 


thickness of the "boundary layer 


P 




pressure 


y 




distance from the wall of the tuhe 


To 




shear stress at the wall 


z 




coordinate in the direction of flow 






The equilibrium condition for the element of the hound- 



ary layer can then he v/ritten as follows (note 16); 

M - pU ^ = - ^ (2r6 - 6^) IT - 2rTr.To (22) 
dz dz dz 

In addition, the following condition appears: Because of 
continuity, the same amount of fluid must pass through all 
cross sections of the tuhe. If, therefore, the layers at 
the wall are retarded, then the velocity of the undi s turlied 
fluid (undisturbed always in the sense that no shear stress 
is transmitted) must increase. Then, according to the equa- 
tion of motion for ideal fluids, this increase in velocity 
must he accompanied by a decrease in pressure, This fur- 
nishes additional boundary equations: 

2 3 

q + U(r - 6 ) 17 = vr TT (23) 



+ pU = 0 (24) 
dz dz 

From equation (2a), u is set into the formi 



and the coefficients a and P are determined from the re- 
quirements : 



u = 0. for y = 0 
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u = U. |a s 0 for y a 6 

■■ " ■■ ; ■ " • •"' ■ oy — ■ -■ ■ - 

Then there results; 



Hesults of Tq from the condition at (25a) must change 
to equation (2) for small values of y: 



•^0 



Consider the relation, which follows from equation (33), 
"between the velocity of the free stream II and the average 
velocity v: 

U = . ^^^^ (27) 

- 32| + 165 

in which ^ = %\ then the steady-state condition yields an 

ordinary differential equation in | , \i?ith variables sepa- 
rated. There is ohtainod then (note (l7): 

.3_100.a 206 . 616^ 14 ^. 
CUI^- 221^ 165 3% '^==(l5^A7i; T^r^^'o (28) 



Instead of determining the quadrature numerically, the fol- 
lowing method of calculation is appli-ed. 

For small values of \ , the higher powers'of ?. can 
he neglected and there is obtained 
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" ■ ""I'f zfr* ■ is' set equal to t, then . | .can 1)6 repre- 
sented as a power series in t, therefore: 

I = At + Bt^ + Ct^ + . . . . , 

Write equation (28) in the foriar 

[1 + p| + qe^]^''* • 

and introduce the above expression; then the fractional 

powers of t drop out and, by developing in powers of t 
and comparing the coefficients on left and right revealing 
the history of the boundary layer at the "beginning of the 
tube, there is obtained: 

i = f = 1-^^ ($) % - (^) (!) * (^) (I) (.9) 

The series is ended at the third term. 

In figure 4 (note 18) | is presented as a function of 

X = ' z ; S = 1 for X =^ 0.686; thus the length of the 

V ^d 

tubeup to the section where the boundary layer fills the 
tube is (note 19): 

z^ = 0.635 d(^^j (30) 

With equation (39) the field of the average velocities 
in the 'ent»rance section is de-termined jLnd.alJ. que for 
which knowledge of its variation is sufficient, can be an- 
swered. Thus, for example, one obtains the resistance of 
the tube between two cross sections at and Zg of the 

initial length; that is, the integral 
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W ~ - 2vn I TQ(z)dz 



"by calculating the difference of momentum transport througlx 
the cross sections and Zg and calculating the pres- 

sure difference times the crossoseotional area and then add- 
ing the two results. 



8. HEAT rHANSFER IN THE ENTEANCE SECTION 



Por the calculation of the temperature field a similar 
consideration is employed by sotting up a heat "balance for 
one element of the "boundary layer. Again, let a warm fluid 
flow into a tuTae with constant v/all temperature ~ 0, 

The only simplifying assumption made is that at a place 
where shear stress is not transmitted, heat transfer will 
not occur. In so doing the small amount of heat which is 
continually carried away at the inner limit of the "boundary 
layer "by molecular condition ig neglected. However, at that 
place the temperature gradient is so small, since the calcu- 
lation is carried out for velocities in excess of the criti- 
cal, that the error committed can "be taken directly into 
consideration. Therefore it is assumed that in the region 
whero the undistur"bod fluid flows with velocity U, the 
temperature alv/ays should "be equal to the entrance tempera- 
ture The heat "balance for the element considered is 

then : 

& 

-i- / uO^ 2Tr (r - y) dy = - 2rTT(i„ (Sl) 

dzJ dz 



Inasmuch as it already has "been seen for hydrodynam- 
ically cODiplete flow that in the case of tur"bulence the tem- 
perature distri"bution is very similar to the velocity dis- 
tri"bution, d is given as J 
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OorreB:|iOnding to the aljove assui^pti ons , the "boundary 
conditions are: 

dab for y » 0 

= « |A = 0 for y = 6 

" oy 

which permit reduction of the three coefficients a, 3, Y 
to a single one. Thus is o"btained 



The first term is identical with the equation for the 
velocity distribution and the second can he interpreted as a 
kind of correction terra to the equation of the velocity dis- 
tribution, This can he easily understood in a physical 
sense. At y ~ 0, temperature and velocity curves begin 
with the same power of y. At y =» 5 , both have horizontal 
tangents. The curves must therefore have a similar charac- 
ter in the intermediate region, 

jPor q.Q it is found that: 

1,340 v^'^^v'^* 1 8/7 + y\ 
<lo " (32a) 



Equation (31) will now furnish, since the variation of 
the boundary-layer thickness 5 is known, an ordinary lin- 
ear differential equation of the first order for determina- 
tion of Y. This differential equation can be brought into 
the form: 

^ + A(z)Y = B(z) 
dz 

The functions A and B are very unwieldy, however, 
so that the general integral of this first order differen- 
tial equation: 
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would reqiiire a very troul)lesome numerical calculation. A 
graphical method is chosenj for it is still pos^i^ble to oTd- 
tain far greater accuracy thsin corresponds to the physical 
assumptions of the problem, If fi/r = | is introduced as a 
new independent variable, the equation is essentially sim- 
plified, the fractional powers of i drop out, and after 
some calculation there is obtained; 

-0.1855 |» + 1.477 - 2.658 | 

y = — — . 

0,1623 I® - 0.701 I® - 23.05 | + 45.4 



-0.369 I* + 3.390 |^ - 16.10 + 37.59 | . , 

— , _ yt (33) 

0.1623 I - 0.701 - 23.05 I + 45.4 



in which Y» = ^ 
<i| 

In this form equation (33) is directly suitable for 
calculation of the directional field of the differential 
eq_uation which is given in figure 5a. The point 1 = 0, 
Y = 0 is the point of origin; all solution curves come from 
plus or minus °o up to a tinique curve which leads to the 
point 0,0, Since Y must likewise be finite for 1-0, 
then for the initial condition there is obtained 

Y = 0 for 1=0 

This is likewise readily understandable from a physical 
viewpoint, .. As long as 6 <g: r, the immediate beginning of 
the tube can differ, either in hydrddynami cal or thermal re- 
spect, from the behavior of a plate in a free stream. It 
will be seen later that the velocity and temperature fields 
are the sane in the case of the plate, Accordingly, at the 
beginning of the tube (z = 0), the temperature distribution 
will coincide with the velocity distributi on. 
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If eq.uation (53) is solved in terms of V and the 
"limiting Vatue tovfard which V • tends- is c-onsidered, for 
i = 0, then there is found: 

lim Y' = -0.032 

|aO . 

All isoclines "begin at the point 0,0; the isocline for 
V = 0,032 with the slope assigned to it, runs into this 
point; however, it is rather weakly concave toward the ah- 
scissa, Tho isoclines of greater slope lie entirely a'bove 
it; those ;*ith a lesser slope lie below it. Prom this "be- 
havior it follows that the solution curve sought must lie 
for its entire length in the narrow strip between the iso- 
cline V = -0,032 and its tangent at the zero point. The 
lino Y = -0.033' I will, thoroforo, represent a first ap- 
proximation with a maxiiaura orror of 12-i- percent. 

In order to obtain a second approximation, Y is set 
equal to -0.032 I + hit) and this expression is introduced 
into equation (33), which then changes into a differential 
equation for h, the family of isoclines of which is shown 
in figure 5"b. In order to increase the accuracy, a thousand- 
fold scale of ordinates is chosen. 

The solution naturally iDegins with h s= 0 for t = 0. 

If e is set equal to 1000 h, then this magnitude, as 
is easily proved hy plotting on a logarithmic scale, is 
given "by the formula e = 1.48 |i-Pe5, so that the following 
expression for Y is finally obtained! 

Y = -0.032 i + 0. 00148 ^'^65 (34) 
This function is given in figure 5a hy the deep, solid 

line. 

In this manner the temporaturo field for the region of 
the simultaneous hydrodynamic and thermal calming length is 
ohtained. Then for d may be written: 

dj^iy^'^li-i. .2.-^1^0^ [0.00148 0,033 nj (35) 

The values ef 6 and | are taken from eqviation (29) 
and figure 4, respectively. For | - 1, ^ "becomes 

^ => Ho (^^^ |l.ll2 " 0.0819 2. „ 0,0305 Cl^^l (35a) 
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that is, the final temperattxre di strihuti on for the hydrody- 
namical s'teady state is not attained- (fig, 6),' which was, 
perhaps, to "be expected 

In order to get further agreement, every solution of 
the differential equation for the temperature field "beyond 
the hydrodynamic calming section, which also satisfies the 
initial condition (ecLUation (35a)), is determined. To this 
end the function represented "by equation (35a) must he de- 
veloped according to the charactetistic functions. (See 
equation (l6),) Since the temperature distribution of equa- 
tion (35a) does not differ very much from the first charac- 
teristic function (equation (16)),. the development is carried 
out with only two memhers, 

The development naturally cannot represent the function 
(35a) quite exactly, "because at one time the velocity dis- 
tri'bution was established with the "influence factor," the 
other time with' a power series development. However, in 
order to o"btain the hest possible transition from the one 
solution to the other, it can "be arranged that the two tem- 
perature curves agree completely in important properties. 

The heat transfer is limited "by the processes at the 
wall; accordingly, it will he stipulated that: (l) both 
curves begin with the same term of the development at the 
wall of the tube, and (2) the flow of energy through the 
whole cross section bo equal. 

This furnishes two equations of condition for determin- 
ing the two coefficients of the development. Hence there is 
obtained for the temperature field in the hydrodynaiaically 
complete region 

^ - jl.OlS e~^^^ [0.9544 x - 0,0213 + 0,0668 x^J 

- 0.051 e'^®'' I- 0.7472 x - 4.275 x^ + 6.022 x^]j> (36) 

Knowledge of the temperature field first shall be used 
for calculation of the unit thermal convective conductance. 
The calming length first will be considered-. 

Equation (36) already has given an expression for, q^. 
The average „j;eraperature of the cross section is then: 

/ *2n(r « y) dy + i>o / 2iT(r - y) dy 
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The expression for cc^^^ (a^ - a, in the hydrodynamioal 
calming region) is then: 



=! 1.340 vC ( -1^ 




1.11^3 - 0.0321 + O.OOlUSl^* 
— ■ i_ TH (37a) 

(H -221+165)^1 |l-0.133l+0.02lH +(0.00ll|S|''^^® -0.032|)(0.520|-0.ll+3l+e"')j 

The continuation . in the second region .' (ag_ = a in 
the hydrodynamical steady-state region) is: 

= 0.03461 vofiLy* 0, 969 e"^^' -f 0, 038 e'^^^ (37^) 
V^/ 0.873 e"^'^'' + 0. 0068 e""^'^ 



The variation of the unit thermal convective conduct- 
ance with location in the calming length of the tube is 
shown in figure 7a (note 20), Write equation (37a) in the 
form 

= KtC (^^^"^^ (kcal/hr '='0 

then the factor K, which is a pure function of i, is 
plotted as ordinate, with i as a^bscissa. Figure 7A (note 
21), in comlJination with figure 4, from which the particular 
values of i can T5e taken, covers all possible cases. In 
figure 7B , with z as the abscissa scale, the variation of 
a for a certain case is shown in comparison with the same 
case for the completely developed hydr odynami-c flow>, It is 
seen that the decrease of a takes place leas ciuickly in 
the first case.. 

9. SUMMAHY AND COMPARISON WITH OB SBRVAT lOK^ 

How comes the point v/here the heat transfer in a tube 
can be surveyed in all particulars. The results will be 
summarized briefly:. 
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With respect to heat transfer, the following cases in 
which the -heat transfer. oteysdiffeE-ent laws are to "be dis- 
tinguished according to the structure of the velocity and 
temporature fields* 

1. Fully developed hydrodynami e and thermal fields »- 

This condition is attained when the fluid has, passed through 
a considerahle portion of the tuhe length, !Phe unit thermal 
convective conductance is constant and is given "by equation 
(31a). 

3, Fully developed hy.drodynamic flow field, temperaturo 

uniform at entranpg .- Heal i zed hy a connected entrance sec- 
tion which is maintained at the original fluid temperature 
hy suitable heating, The unit thermal convective conduct- 
ance is dependent upon the location in the tube, falls very 
quickly from its maximum value, and asymptotically approaches 
a constant minimum value, (See equation (21),) 

g. Uniform velocity and temperature di strihuti ons 
across the section at entrance ,-* The unit thermal convective 
conductance is likewise dependent upon the location in the 
tuhe, hut falls to a minimum value more slowly. The point 
unit thermal convective conductance is given hy equations 
(37a and 37h ) , 

4, The application of heat "begins at a section somewhoro 
in tho middle of the calming length ,- For this last case a 
good approximation for the unit thermal convective condxict- 
ance a is obtained "by drawing the curves which represent 
the variation of a with the location in the tube for cases . 
3 and S, in the same system of coordinates hut with the saero 
point of the ahsoissa scale for (a fully developed hy- 

drodynamic flow) displaced hy the distance I "between the 
"beginning of heating and the inlet section. Since the curve 
for has a much steeper slope, it will cut the curve for 

Ot^; the envelope (note 32) represents (to a first approxima- 
tion) the a distribution for this special case. 

The great differences in the results of the individual 
experimental works are now understandable, Whereas Nusselt 
ascribed this, at the conclusion of his \irork on heat trans- 
fer in laminar flow, exclusively to the conditions mentioned 
under point 2, now the possibility of a series of factors 
which influence- the process by interohangeabj-e combinations 
may be seen. 
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The two most careful investigations known, those of 
ITusselt and Jordan (references 6 and 10, respectively), used 
a right-angle gas approach to the measuring section. At the 
"beginning of this section the flow was not completely devel-^ 
oped. OTjviously case 4 is to "he considered here. Ho fur- 
ther experimental data, however, havo appeared to dato. The 
■beginning of the thermal action, as woll as the exact posi- 
tion of tho first tomporature-measuring station, cannot ho 
ascertained accurately,* An exact evaluation of the exper- 
imental results on the hasis of the a"bove-mentioned theory 
is not possihlo for these oxpotiments, Hovortholoss , a 
series of experiments "by Nusselt were investigated to deter- 
mine the magnitude of the measured unit thermal conductances 
with respect to tho minimum a. The results arc compiled in 
tahle 2, Considering the dimensions of Nusselt's apparatus, 
it is seen that these figures are affected hy them, which is 
to he expected according to tho ahovo-mentionod derivations, 

TABLE 3 



Exper- 
iment 
nuQ- 
ber 


I 


Pm 


1 ' ' ".' 
^m 




V 


•^eas 


Omin 


Differ- 
ence 

(percent) 


7 


0.5906 


1.153 


^3.7 


1.2lf5 


S.52 


35-7 


33.1*3 


e,h 


10 


.59S3 


1.152 


35.3 


1.27s 


IS. 33 


65.5 


60. U 




13 


• 5880 


1.171 


31.6 


1.315 


27.2 


91.7 


82.6 


11.0 


19 


.7S7O 


2.050 


69.0 


2.050 




57.3 


50.9 


12.5 


2k 


.6022 


I.SS3 


30. U 




29.9^ 


1H6.6 


127.2 


15.3 


30 


.6178 


3.96 


31.9 


1+.H5 


10. S7 


12li.S 


103.8 


20.0 


In 




6.97 


29.2 


7.91 


9.27 


162.9 


iin.6 


15.1 


5^ 


.5S63 


9.9s 


26.0 


11.1+1 


lO.OU 


233.0 


19S.0 


17.7 



If experiment 10, for example, wore evaluated with tho 
assumption that case 3 is ^0 he considered hero, which comes 
Closo to tho experimental conditions in any case, then a 

♦These results are hased on a letter sent to Latzko hy 
Prof, Husselt. 
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value of a is o^btained, which is 10,2 percent a'bove a^^^ 



tion of velocity at the first tomperat'p.re'-measuriiig station 
was not complotoly fulfilled, so that tho value of a must 
fall sonexirhat lower; tho moasurod value is 8,4 percent above 

*^min» 

According to equation (29^), the length of the calming 
section is proportional to ^/vp', accordingly the ratio of 
tlie calming length to the total length of the measuring sta- 
tion must increase with increasing velocity, and density and 
therefore tho values of a also must increase somewhat. 
Experiments 7, 10, and 13, in which the velocity is varied, 
the other parameters being maintained reasonably constant, 
show this clearly, as do experiments 7, 19, 41, and 54, in 
which the values of p are changed up to a ratio of 1:10, 



Finally the q.uestion of the practical computation of 
the amount of heat transferred in tubes must be discussed. 
Because of the fact that tv/o different solutions for the 
temperature field have been obtained, depending upon v/hethcr 
or not the hydrodynamic field is fully developed, it is aec- 
essarj?" also, in the computation of the amount of heat trans- 
ferred in a certain section of the tube, to determine in 
which region the flow takes place. However, several gener- 
al remarks, which hold for both regions, must be presented 
first. 

The amount of heat transferred between two cross sec- 
tions at an interval I » - Zg can always be found in 
two ways, 

1, The calculation is referred to the volume of fluid 
passing through the cross section per unit time. Since the 
velocity and temperature distributions are known, then tho 
flow of heat which passes through tho cross-sectional area 
per unit of time can be obtained by integration. The amount 
of heat transferred at a strip of wall of length I between 
Zx and Zg thus is given by the difference of two integrals 




10, PRACTICAL COMPUTATION 0¥ THE AMOUNT 



OF HEAT TRANSFBRRBD 



q= Oy u(y,zj)*(y,zj)df - C ^ u(y,za)d(y,za)clf 

f ■ ' f 




(38) 
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where df is an element of area in the cross section- the 
-f-oipm-o-f— which depends upon the choice of y. 

2i Only the processes at the wall are considered. I'or 
this purpose the equations for q,o liave heen set up. With 

the foregoing notation, the amount of heat transferred is 
then: 

Q = 2rtr ^(i^(z)dz (39) 
za 

the term in qo replaced Tjy the difference at the 

cross section hotween the wall temperature and the 

temperature at the axis of the tuhe. 

The expressions for q^ are, however, rather unwieldy, 

so that equation (38) is usually preferahle. As an example, 
equation (38) will he applied to the most important cases. 

la. Por the case where a uniform temperature distrihu- 
tion and a uniform velocity di strihuti on prevail in the ini- 
tial cross section, the heat transfer in the entire entrance 
section is: 

Q s; 0.115 r^TTvOdo (40) 

The Value is constant since the temperature and velocity 
distrihutions coincide at the initial and the final cross 
sections. The special cases differ in the length of the 
tuho section from the inlot opening to the fully developed 
hydrodynamic state. 

Ih, Considiering a section of the tube from the inlet to 
a cross section at a distance L and letting I he the 
length of the entrance section, then, (using equation (36)), 
the total heat transferred in the length Ir is: 

Q =* r\v-Odo-|^l - 0.886 e"^^-^^- - 0.0037 e'^^^^j^ (40a) 

whore li - L - I, 

2, In the fully developed hydrodynamic state the amount 
of heat transferred from the beginning of the heating (z=0) 
to a cross section z at a distance I from it, comes to: 
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Qr: Trr^vO^o - 0.985 e~^i ^-0.013 e"^2 1^0, 0022 e'^^^^^ (41) 

11. BSTIMAT.IOIT OP THE IITTIRSTAL HEATING 
or THE PLUID BY fHIOTIOH 



In the previous derivations the generation of heat in 
the interior of fluids "by friction was entirely disregarded. 
An exact consideration of this is not possible as long as 
the pulsating velocities are not known individually, knowl- 
edge of the distribution of the (time) average velocity be- 
ing espocially insufficient for this, since the pulsation 
velocities contain the dissipation function in quadrature 
terms. However, an attempt will be made to determine, by an 
approximate consideration, in which velocity region the 
above-mentioned neglect of this term is allowable. In doing 
this, only processes in the fully developed hydrodynamic and 
thermal states- are considered; that is, the unit thermal con- 
ductance must be independent of the location in the tubo. 

It is assumed, for example, that a cold fluid flows 
through a heated tube, an element of volume bounded by the 
tube and txiro cross sections at a distance dz is considered, 
let be the constant wall temperature; let d(z) be the 

average temperature of the cross section at the position z; 
and let the difference between vrall and average temperature 
in the initial cross section be designated by ^ q. Then the 

heat balance on the element reads! 

r^TTvOd* = a(*^ - i&)3rTrdz + T t,r^ 2vnd.z (42) 

where T is the thermal equivalent of work and t, is the 
ratio of the frictional resistance to the square of the ve- 
locity* Equation (42) then is written in the form: 



— + & = ( d„ + T 

dz rvC 




-A 

a / 



1 1 



The term J ii- has dimensions of temperature; 

a 



V will be called the friction temperature and will be 
a 



designated by ^jj^. 
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"Wit-h -ao-as-ideratton, o^;,t^^^ at z = 0, 

the solution of the differential equatloh (43) is: 



Kence, the temperature of the fluid is proportional to 
the wall temperature increased "by the friction temperature. 

If this relation is plotted in a system of coordinates 
with 9 as ordinate and z as abscissa, then there is. oh- 
tained the clear results of figure 8, which are drawn for 
the case of air. flov/ing with 

V s 300 m/sec r = 0,035 m p = 1 atm 

0 = 0.383 kcal/m^ v = 0.175 cm^/geo 



The fluid temperature asymptotically approaches a limit 
which is eqiial to the ^vall temperature plus the friction 
temperature. Since the aniount of heat transferred is pro- 
portional to the areas (xirhich are crosshatched in fig. 8) 
"betxtfeen the straight line = constant and the curve of 

the temperature of the fluid, it is seen that there exists 
such a relation "between the length of the tuhe section and 
the velocity of the stream that, for a given length, there 
exists a certain velocity and, for a given velocity, there 
exists a certain length of tuhe for which the maximum heat 
is transferred, 

t 

The factor is then? 

a. 

^ = 0,103 
a vc 

and, therefore, if the slight variation of the heat capac- 
ity with temperature is disregarded, it is a function of ve- 
J^ocity alone. for c = 0.338 kilocalorie per kilogram °0, 
for air has Vepri calculated fC'r-e'e^^ral v el. 0 cities and 
compiled in the following tahlo: 

V = 10 25 50 100 150 300 m/sec 

ft 

*H 0.103^ 0.634° 3.54"^ 10.15° 33.80° 40.6° 0 
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i-t can "be seen that neglecting ■6^ is directly permis- 
siTjle in most practical cases. 

2he heat carried away is obtained by .substitution of 
(43) into the equation: 

dQ s aSrTTdz - ■dg) 

and integration between 0 and Zs = I to 

q. = ^or^TTvC \^1 - e rvTTy « ^ ftj^r^rrvG (44) 

in which higher powers of are neglected. 

12. HEAT TRAHSI'ER OK A PLAT PLATH 



As a second geometric configuration for which the heat 
transfer will be calculated, the flat plate parallel to the 
direction of flow is chosen. The plate shall be so thin 
that the influence of the forward edge can "be neglected. 
Prandtl, in the previously mentioned work, has already shown 
that for the case of an infinite thin plate, v/hich is moved 
parallel to itself through a fluid, the velocity field and 
temperature field agree, if the hoat from the internal fric- 
tion is neglected. 

If u denotes the velocity vector and p denotes the 
pressure, then, for the time change of the momentum vector 
of an incompressible fluid referred to a unit of voluno, the 
result is, neglecting gravity effects (note 23): 



Du 

P jT- = ~ grad p + nAu 



dt 

I) 
dt 



(45) 



bt 



u -2_ + 



On the other hand, the Kirohhoff differential eoLuation for 
the temperature field is: 
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■ ^ ' - — ■ ■■ -^^Si ^ Ji. Ad (46) 

dt . cp ■ ■" — - ■ 

ffor the plate, p is a constant, in case it is infi- 
nitely thin; or, in case the plate thickness is finite, the 
average value of p over a certain region is still a con- 
staiit; grad ' p is, therefore, equal to 0 and it may "be soon 
that the two equations (45) and (46) agree in form. When 

~ = 1, a solution of (45) also will "be a solution of (46), 

C|J> 

If use is made of this condition when the- solution for 
the velocity field, as given in the previous work Tjy 
Ton Karman, is accepted, thero results i 

u = U 8 = 0,366 (x^)''^ z (47) 



and correspondingly: 



The condition for thermal equili'brium in an element of 
the "boundary layer is J 



dz 



<lo 



Introducing equation (47), solve for q^ and ohtain: 

* q^ = 0.0285 *o OU (^^""^^ (48) 

The total amount of heat leaving a plate strip of unit 
width is then (note 24) ; 




(I = length in the direction of flow). 
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13, !rEAKSMISSI01ir lU LICiUIDS 



In the preceding section it has "been seen that the laws 
for the heat transfer in a turbulent fluid stream, derived 
for the two most important Tjasic geometric forms, lead, in 
the case of flowing gasos (and superheated steam), to re- 
sults which also agree well q^uantitatively with experiments. 
Here the statistical hasis of the kinetic theory of gases 
can he used to arrive at a uniform concept of the molecular 
prooossos as well as the eddy prpcossoe in friction, on the 
one hand, and in heat transfer on the other hand, 

Naturally, one cannot transfer the foregoing simple 
considerations directly to liquids, where the effect of the 
molecular forces of cohesion may no longer he neglected, 
Whoreas the heat and momentum convection through the eddy 
system also represent here processes which are similar in 
character, this is no longer true of the molecular conduc- 
tion of heat and momentum, which finds its expression in. 

Xo 

that the ratio is very different from 1, Per water, 

C|j. 

the magnitude of this ratio, which is quite dependent upon 
temperature, is ah out 0,1, 

The mutual law for the molecular and eddy phenomena of 
internal friction, which is represented hy tho coefficient 
of turhulence, will ho applicahlo to tho propagation of heat 
only in that region in which eddy processes dominate. How- 
ever, this is the case for the entire mass of fluid up to a 
very thin layer* at the wall. Accordingly, the differential 
equations derived in the foregoing will maintain their va- 
lidity everywhere except in this thin layer. 

To attain a suitable description of the heat transfer 
in fluids, it will he necessary to seek a transition from 



♦Closer investigation shows that this layer is much 
smaller than the boundary layer itself, which was defined as 
the region in which shear stresses are transmitted, Tor the 
right circular tube, the thickness of this layer in the hy- 
drodynamricaily perfect state is given by^the ejcpression 

6 « 6,51 -«r- where d is the diameter and H is the 
B,eynolds number. 
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the statistical law for the interior of the mass of fluid 
(coefficient of turljulence ) to the molecular law in the .very 
neigh'borhood of the wall (thermal conductivity). Mathemati- 
cally, . this transition can "be made hy a modification of the 
"boundary conditions. This extension of tho theory shortly 
will "be discussed moro closely. Still, it is to ho noted 
that oven the previous experimental results of research, 
which relate exclusively to heat transfer in water, show re- 
sults deviating so auoJx from one another that only with dif- 
ficulty can a picture of the process "be made clear to some 
degree. 

At considerahle expenditure, Soenneoken (l91l) under- 
took experiments on the heat transfer of virater in tuhes, 
These experiments are freq^uontXy oitod in modern literature 
as authoritative. His results are summarized in two formu- . 
las for the unit thermal conductance a-J 

1, Smooth surfaces: 

a = 2020 (1 + 0.014 Ti ) ■■ ■ ^^.g^-v -^ 

d • hr ffl^ °C 

S, Rough surfaces: 

v°*'' kcal 
a = 735 -573- (1 + 0.014 T., ) 



d " hr m^ °0 

where 

V water velocity, meters per second 
d tuhe diameter, meters 

2?j^ internal tuhe-wall temperature, degrees centigrade 

By "smooth surface" is understood a seamless drawn-hrass 
tuhe; whereas the experiments with "rough surfaces" were 
performed with iron tuhes. These formulas directly contra- 
dict the fundamental ideas on tho nature of the heat transfer 
presented hero. It is out of the question that the unit 
thermal conductance can he smaller in the case of rough sur- 
faces than it is for smooth surfaces. Because of the in- 
creased eddy formation, the eddy transport is increased and 
conseauently the unit thermal conductance is increased, Tho 
smaller values which were measured in tho case of tho iron 
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tubes are , most prolDably ascribalile to an incorrect determina- 
tion of the actual wall temperature due to. tJie presence, of 
layers of Tjoiler scale and rust. 

Extensive experiments on condensers were performed 'oy 
Josse in OharlottenTjurg, Josse measured the over-all unit 
thermal conductance k of condensing steam to w«^tor, which 
is defined "by the formula (reference 11): 




where 



ai,aa unit thermal conductances of the fluids 
X thermal conductivity of the partition 

6 thickness of the partition 

Since Josse substitutes for X and 6 the known values and 
uses the figures of ITichol (about 4500 kcal/hr °0) , which 
represent in magnitude a raean between the figures of 
Soenneoken for smooth and for rough surfaces, for the heat 
transfer in tho water, he obtains unusually high values for 
the heat transfer from condensing steam to metals. Accord- 
ing to these experiments the unit thermal conductance for 
condensing steam would be about seven times as large as that 
for flowing water. 

These results are likewise not understandable. Accord- 
ing to Husselt (reference 12), the process of condensation 
on a cold perpendicular wall produces on the cold side of 
the surface a film of water in which occurs all the drop 
from stoam to \taLl temperature. The film of water clings to 
tho Wall, the remaining layers floy downward under the ef- 
fect of gravity. Since the thickness of the film of water 
is very small in any case - fractions of a millimeter accord- 
ing to JTusselt's calculation - the flow obviously must be 
laminar. It is impossible to understand why the unit thermal 
conductasice from water to th© metallic, wall should be seven 
times greater on tho ono side of tho wall (whorb tho state 
of flow is laminar) than on tho other side (whoro there is 
turbulent flov/). 



HACA TU Ho, 1068 



40 



.Ih.,e. expjBrime^^^ Josse accordingly ought to 

have soae other explanation so that k, the over-all unit 
thermal conductance, is represented "by two terms, approxi- 
mately equal in magnitude, -i- and -i^ where aj, and ag 
are average vailues. 

It seems quite prohahle that, for this reason, the unit 

thermal conductance for water has higher values than were 
frequently assumed previously. An extensive experimental 
investigation of the heat transfer to fluids seems to "be 
urgently needed in order to bo able to test tho accuracy of 
the thoorotical calculations. The difficulties which are 
encountered due to formation of rust and scale when water is 
used, suggest the use of other fluids, such as oil, for ex- 
ample* 



Translation Tjy 1. M. K. Boelter, 
<3r. Young, and A. 0, Guibert. 
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- EXPLAITATORy .ITQfBS BY f HE TEANSLATORS 



NOTE 1 



The eq^uation, 



Can "be obtained in the following manner. Postulate that the 
l/7--po\irer equation for the velocity distriliution holds near 
the wall and that the shear stress at, and in the vicinity 
of, the Wall is a function of y, the distance from the wall, 

and of the velocity gradient at that point, Then, 



divide hy the particular velocity gradient, which gives 




T 




6y 



Prom the l/7~power equation (see note S), 




there is olitained: 




or 
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In the -vicinity near -the wall, let 

'° (-If) =--0 (-If) 

therefore 



or 




Latzko dropped the subscript (y-*.o) and stipulated that the 
equation is valid only in the vicinity of the wall. 



NOTE 2 



Sydney Goldstein (Modern Developments in Fluid Dynamics, 
Clarendon Press (Oxford), 1938, pp. 339-340) reoommends a 
value of B = 8,7, 



NOTE 3 



Por fluid flow near the wall in smooth pipes the velo- 
city u is determined by the following variables: T^, the 

shear stress at the wall; y, the distance from the wall; 
p, the density of the fluid; and u, the kinematic viscos- 
ity of the fluid. By dimensional analysis there is obtained 




A comparison of this equation with the Blasius resistance 
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formula (empirical) for flow in smooth pipes (see Goldstein, 
Modern Development in iFluid Dynaini 339-340) allows 

the determination of the magnitude of the exponent n and 
K; so there results the l/7-power equation for velocity 
distrihution. 




HOSE 4 



As an attempt to correlate the hydrodynamic principles 

presented "by latzko with the more recent knowledge of veloc- 
ity and shear distributions, the Prandtl mixing length was 
calculated from equation (?) and compared with those derived 
from Prandtl 's and Zarman ' s logarithmic formulas for veloc- 
ity distribution and also with that obtained from Hikuradse's 
data.* Although there are some inconsistencies in latzko's 
equation (inconsistencies also appear in Prandtl 's and 
Karman's equations), the variation of the mixing length does 
have the same trend as that calculated from Nikuradse's ex- 
perimental data and differs from it only by a constant ratio 

of approximately 1,25 Q^i^^^^f^'^f-^ - 1.25^. The method of 

Calculation and the results are shown in table A-1 and fig- 
ure A-1, 

The most important inconsistency appears in the deter- 
mination of the velocity distribution near the wall. Recent 
developments indicate that the velocity is a linear function 

"^o ^ 

= ^ y J of the distance from the wall in the laminar 

sublayer, Latzko<s expression for the velocity distribution, 
however, approaches the l/7-power equation near the wall and 
would yield an infinite instead of a finite velocity gradi- 



ent at the wall 




I#atzko..,may be. justified in 



using such an expression by the fact that his expression for 



♦Seo Bakhmetev, Boris A, ? The Mechanics of Turbulent 
Plow. Princeton Univ, Press, 1941, p. 73. 
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shear stress is ir = K (y) which could yield a finite 

oy 

value, for Tq, however, ■because the product of ziero and in- 
finity is indeterminate. 



The term Pwx , called hers the coefficient of turTju- 
lence, is identical with the modern term "eddy dif fusivity 

e in ■ 



From this point on, Latzko uses y for the distance 
from the pipe axis. Up to this point, however, he used y 
as the distance from the wall and y for the distance from 
the axis. To avoid confusion, this translation continues to 
use y for the distance from the wall and y for the dis- 
tance from tho axis only. Tho equations that follow havo 
Tieen altered (from the original) to conform with these orig- 
inal definitions. 



NOTE 5 



T a pe 



NOTE 6 



NOTE 7 



The original article gave the equation, 




It should he 




NOTE 8 



The derivation of 
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- ............. ^4, / a _Sv% 5.^ 

^ (2r)%3 V Sr / ay 

from equations (8), (2) , (3), and <4). 



Equation (8): 



Equation (2) : , 



7 

B 



Equation (s): 



Equation (4): 



u = u 



max 



-CD 



Erom equation (3) there is oTjtained 



u 



u 



B 



But by equations (3) and (4), 



or 



u 



therefore 
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SuTJstltuting into equation (8), 



(0 



^ W 



and since 



\ 2r / by I y 



(7X8) - • - - ^- - ^ 



1 - ( 

y ~ r^y 2r 



thus 



(7 X 8 ) 4 V V 0 / r " y N ^ ( 1 - X- ) 
b'^ (2r)^a8 V 2r / by ^ 2r / 



and taking ( 1 - f^j « 1, which is approximately true 
the vicinity of the wall (y small), yields 



0.199 V % 0 / r^ ~ y^N 
' (2r)^28 V 2r y 



by 



EOTB 9 



transformation of equation (ll) to equation (12). 
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BqLUation (ll): 



dz 




where 



Let 



then 



„ V4 . ,%8 

8 V (3r) 
7 X 0,199 V 



/_y -kz 
= g (y) e 



and 



kz dg 

- S3 e 
&y dy 



|1 = _ ke-^^g(sr) 
0 z 



Substituting in equation (ll) yields 



y** = (1 - x^) 



To simplify, let 



or 



then 



2y dy = - 7 x®dx 



dy = - 



2y 

B;r suhstituting, there is obtained 
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where 



w = 49 kK 



(I) 



HOT25 10 
The original article gave 



8, 

w = 4gkK 



It should read 



w = 49kK 



rroTB 11 



(a) Sorae references to the Ritz method of solution of 
differential equations. 

Sitz, Walter: UTser eine neue Methode zur LOsung 
gewisser Variationsprobleme der mathematisohen 
Physik. J. f. reine u. angew. Math.., vol. 135, 
1908 , pp. 1-61. . 

Uher eine netio Methode zur LBsung gewiisser Rand- 
tvertaufgaheh. GOttinger Hachrioh.ten , 1908, pp. 
236-248, 

Both in. "Oeuvres de Walther Ritz," Paris, Gauthier- 
Villars, 1911. 

ilimoshenko, S. : Vihration Prohlems in Engineering, 
D. Van Uostrand Co., Hew York, 193&, p. 259. 

(h) The manipulations involved in the Ritz method, as ap- 
plied to this proTjlenii 

• 1, Insert the Logendre spherical functions of the 
first kind into equation (14). Thon, 
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3. Now differentiate equation (14) with respect to x 
to obtain ^ 

3. Substitute into equation (is) , the value of g and 

r— obtained in steps (l) and (2) and integrate, 
ox 

1 Then, 



+ 



J \8 loA U 1+ g (S)2 J 

/iH w/s 15 3^5.15 70 210 70.15 



2 



h. Then differentiate the eqijation obtained in step (3) with 

respect to gi, gs, and g,, yieldisag the three equations 

aj 8J aj 

for r — ' r — ' and r — where 



and set them equal to ^ero. 



1^ =0= (1.75-0.2 w)gi+ (O.S75r. 0,1167 w)g2+ (1.22- O.Ol^ w)g3 



UAOA TM Ho. 1068 51 



= 0= (O.S75- 0.1167 w) ea+ (6.5^-0.092 w) ea+ (3.i^3-o.o6 w) g. 



a J 

§^=;0= (0.22-0.0U2 w) ei+ (2.1+1-0.05 w) g3+(u.U5~ 0.053 w) 

5. Set = = = 0; Tiy setting the -determi- 

oSi ogs 

nants of the three eq.-uations to zero, a third 
degroo algctiraic eq.uation involving w is olj- 
tainod; 

109,3 - 13.3 w + 0.088 - 0. 0000S7 w'* = 0 

6, Solve this third-degree equation in w, to ohtain 

the three characteristic values, w^^ = 8,713, 
wg = 164.36, and W3 = 1700,40 as given hy 
Lat zko . 

7» Substitute these characteristic values into the 

three equations in step (4) to obtain three sets 
of Si » Spt and g^ ; and insert these into 
ecLuation (l4) to obtain ecLuation (I6) and thus 
equation (I8). 

8, The values of ai , aa » and aj in equation (I8) 
are obtained by the method of least squares, 
thus yielding equation (l9). 



NOTE 12 

References to Calculus of Variation: 

1, Woods, Frederick S. : Advanced Calculus, Ginn & Co., 

19S6. 

2. Bliss, Gilbert A.: Calculus of Variations, 1925. 

IT 0TB IS 



The function 
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satisfies 25uler's differential equation, which is the neces- 
sary though not sufficient condition for a minimum, 



6f _ A- 

&g dx \dgV 



for 



and 



^ = - 2wx^g 
og 



dx \&g»/ dx I. \dx/J 



therefore 



il - ^ . Iswx'g - 2 -i- ((1 - x') 

bg dx \6gV dx I Vdx, 




and since 



A. Ui > ,^7) ^] = -wx'g (12) 
dx L dx J 



thus. 



5f ^ , 0 



&g hx \hg 
NOTE 14 



References to Legendre • s Polynomials. 

1» Jahnke, , and Emde, P.; Tunktionentaf eln mit formeln, 
und kurven, Dover PulJlications , 1943, 



2. 



Woods, S, ; Advanced Calculus. Ginn & Co. , 1926, 



ISTACA m TSo, 1068 



5S 



NOIE 15 



K. ten Bosch in "Die Warme-ttTjertragTing" (Julius Springer 
(Berlin), 1936) rearranged equation (21) to yield: 



-a , 7 



TP r, ^ 0.1 



•••)■ 



which when more terms are added "becomes 

" ~ 3 « £ -'39.37 -31.96 

0.03SU vC J ^ 'oTisd o.a5 d 0.35 d 

= ^ 6.35" Si + 0.1 c"® + 0.9 e - 0.023 e > 



a= 5,36.35 

which is a more convenient form of equation (21). Notice, 
however, that these equations do not yield an infinite unit 
thermal conductance at the entrance v;here z - 0, which is 
in contradiction to Latako's statement in the paragraph fol- 
lowing equation (2l), 

NOTE 16 

The momentum equation. of the boundary layer can he 
written 

Referring to figure 3, it can he determined that: 

6 



^^^''(h /p-^^rr (r-y) dy) 



dz 



0 

the flux of momentum across cd that exceeds that across 

ah 

dJ 
dz 



= ^ nv^'o fu^ / 1374 „ 2503 {\ 
^ TTr p ^207 345 

165v (8| - 22) ri274 , 2303 

2u (_ I i \ 

(4| - 221+ 165) ^207 690 



d| bq. 
-i- - pU — dz 

dz dz 
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the -im^ard flux- of mojaen turn across . "bd 

- pU 1^ =+ 2pUnr^|^U (1-1) 

165v (Bj^ 22)(1- 31+ i^) 1 dE dp . 3, 

" 7-75 — (2rS - 5 ) dz 

(4|** - 32| + 165)^ J dz dz 

difference in pressure "between aB" and cS" 

- ^ (2r6 - 6") = . 2^rapTJ(3^-^")165y(4^,> 11) ^L^^^^r 

(4P ^ 22| + 165)^ dz ° 

retarding force at th& wall 



where 8 

Q = j/^ 2Tr (r - y) u dy 

o 

By substituting these values into etiuation (22), there is 
obtained equation (28), 

ITOa?! 17 

Though the original article gave z in equation (28), 
it is obviously Zq after considering the limits of the 

integration, 

JSrOTB 18 

The original article indicates that the curve in figure 
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4 is a graph of equation (29), However, on replotting the 
expression, a discrepancy was found. It appears (see fig, 
A-fll) that the curve given "by Latzko in figure 4 is in er- 
ror, 

NOTE 19 

Srom figure A-II, it is seen that for 1=1, X = 0,686 
.vd/ \d 

V 

or 5/. y^jiS. '4. 



= 0.686 



Zo = (0.686) d f— j 
= 0.635 d (^I^y* 

Latzko derived from figure 4, that 
Zo = 0,693 Cl {^^* 



HOTS 20 




V 



This plot is given in the original article under fig- 
ure 7A with ordinate misrepresontod as K. 



NOTE 21 



A plot of K against I, 
■by figure 7A, is missing in the 
presented here in figure A-III, 



which was to have "been given 
original article. It is 
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,. N.OI.B ,22.. 



An example of case 4. 



^B-a for case 2 

cL 



ct^ = a for case 3 
= a for C8,se 4 




The original article gave the ecLuation, 



It should "be 



Du 

p ss - grad p + \jAm 
ds 



Du 

^ dt" " " ^^^^ P 
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KOTB 24 



Kxis equation is derived for a fluid the Prandtl number 
of wliich is unity, GollDurn* gives an expression for heat 
transfer at plane surfaces which is valid for Pr other 
than unity, 

Oolljum's equation is 



which is practically the same as Latzko's equation for Pr 

=: 1 



*Con3urn, Allan : A Method of Correlating Forced 
ConvGction Heat transfer Data and a Comparison with Pluid 
Friction, Trans, As. Inst. Chejn, Eng., vol, XXIX, 1933, p. 
139, equation (22). 
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Figure lA, 



Pig. 1 



Figure lA., The three 
character- 
istic! functions a, b, 
c plotted against x as 
abscissa. 



Pigare IB.- The same 

plotted 
against y as abscissa, 
d = initial distribu- 
tion for z' = 0. 




Figure IB. 
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?igure 2,- Variation of heat transfer coefficient along pipe length for 

hydro iynamically fully developed state (c = 0.304 Cal/m3, v = 
18.3 m/sec, v= 0.175 cm2/sec, d = 2.2 cm). 




figure 4. 
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Figs. 6, 
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Figure 6.- End temperature distribution for hydroclynamically developed 

state (e) and distriliution at the end cross-section of the 
entrance run (h) « 
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Figure 7A. 
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Figure 7B. 
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